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Abstr act. When a plane shock hits a wedge head on, it experiences a reflection- 
diffraction process, and then a self-similar reflected shock moves outward as the 
original shock moves forward in time. The complexity of reflection-diffraction 
configurations was first reported by Ernst Mach in 1878, and experimental, 
computational, and asymptotic analysis has shown that various patterns of 
shock reflection-diffraction configurations may occur, including regular reflec- 
tion and Mach reflection. In this paper we start with various shock reflection- 
diffraction phenomena, their fundamental scientific issues, and their theoretical 
roles as building blocks and asymptotic attractors of general solutions in the 
mathematical theory of multidimensional hyperbolic systems of conservation 
laws. Then we describe how the global problem of shock reflection-diffraction 
by a wedge can be formulated as a free boundary problem for nonlinear con- 
servation laws of mixed-composite hyperbolic-elliptic type. Finally we discuss 
some recent developments in attacking the shock reflection-diffraction problem, 
including the existence, stability, and regularity of global regular reflection- 
diffraction solutions. The approach includes techniques to handle free bound- 
ary problems, degenerate elliptic equations, and corner singularities, which 
is highly motivated by experimental, computational, and asymptotic results. 
Further trends and open problems in this direction are also addressed. 



1. Introduction 

Shock waves occur in many physical situations in nature. For example, shock 
waves can be produced by solar winds (bow shocks), supersonic or near sonic air- 
crafts (transonic shocks around the body), explosions (blast waves), and various 
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natural processes. When such a shock hits an obstacle (steady or flying), shock 
reflection-diffraction phenomena occur. One of the most important problems in 
mathematical fluid dynamics is the problem of shock reflection-diffraction by a 
wedge. When the plane shock hits a wedge head on, it experiences a reflection- 
diffraction process, and then a fundamental question is what types of wave patterns 
of reflection-diffraction configurations it may form around the wedge. 

The complexity of reflection-diffraction configurations was first reported by 
Ernst Mach [109] in 1878, who first observed two patterns of reflection-diffraction 
configurations: regular reflection (two-shock configuration) and Mach reflection 
(three-shock configuration); also see [5j 1961 1118) . The problem remained dormant 
until the 1940's when von Neumann, Friedrichs, Bethe, as well as many experi- 
mental scientists, among others, began extensive research into all aspects of shock 
reflection-diffraction phenomena, due to its importance in applications. See von 
Neumann [142l [135] and Ben-Dor [5J ; also see IM HI H [95J H051 [T5TI IT52] 
and the references cited therein. It has been found that the situation is much 
more complicated than what Mach originally observed: The Mach reflection can 
be further divided into more specific sub-patterns, and various other patterns of 
shock reflection-diffraction may occur such as the von Neumann reflection and the 
Guderley reflection; see [5j [M [77J USB US IT291 [T55l IT561 [1571 IT591 [T42l IT45] 
and the references cited therein. 

Then the fundamental scientific issues include: 

(i) Structure of the shock reflection- diffraction configurations; 

(ii) Transition criteria between the different patterns of shock reflection- diffraction 
configurations; 

(iii) Dependence of the patterns upon the physical parameters such as the wedge 
angle 9 W , the incident- shock-wave Mach number M s , and the adiabatic exponent 
7 > 1. 

Careful asymptotic analysis has been made for various reflection-diffraction 
configurations in Lighthill [104] . Keller-Blank [93] , Hunter-Keller [90] . Morawetz 
[TT6] . [671 [T54l l84l [871 IT421 IT45] . and the references cited therein; also see 
Glimm-Majda [77 . Large or small scale numerical simulations have been also 
made; see, e.g. [5j[77], |SB EH EH S3 H2DI IDS] , [51 1521 ITU 1721 ITS! 1521 [144] . 
and the references cited therein. 

However, most of the fundamental issues for shock reflection-diffraction phe- 
nomena have not been understood, especially the global structure and transition 
of different patterns of shock reflection-diffraction configurations. This is partially 
because physical and numerical experiments are hampered by various difficulties 
and have not been able to select the correct transition criteria between different 
patterns. In particular, numerical dissipation or physical viscosity smear the shocks 
and cause boundary layers that interact with the reflection-diffraction patterns and 
can cause spurious Mach steams; cf. Woodward- Colella [144] . Furthermore, some 
difference between two different patterns are only fractions of a degree apart (e.g., 
see Fig. 5 below), a resolution even by sophisticated modern experiments (e.g. 
107J) has been unable to reach, as pointed out by Ben-Dor in [5J: 11 For this reason 
it is almost impossible to distinguish experimentally between the sonic and detach- 
ment criteria!'' (cf. Section 5 below). In this regard, it seems that the ideal approach 
to understand fully the shock reflection-diffraction phenomena, especially the tran- 
sition criteria, is still via rigorous mathematical analysis. To achieve this, it is 
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essential to establish first the global existence, regularity, and structural stability 
of solutions of the shock reflection-diffraction problem. 

On the other hand, shock reflection-diffraction configurations are the core con- 
figurations in the structure of global solutions of the two-dimensional Riemann 
problem for hyperbolic conservation laws; while the Riemann solutions are build- 
ing blocks and local structure of general solutions and determine global attractors 
and asymptotic states of entropy solutions, as time tends to infinity, for multi- 
dimensional hyperbolic systems of conservation laws. See [181 I19L 1201 1211 1771 
WE 13 [991 rrUTl [mi fT24l [T23l IT56] and the references cited therein. In this 
sense, we have to understand the shock reflection-diffraction phenomena, in order 
to understand fully entropy solutions to multidimensional hyperbolic systems of 
conservation laws. 

In this paper, we first formulate the shock reflection-diffraction problem into 
an initial-boundary value problem in Section 2. Then we employ the essential 
feature of self-similarity of the initial-boundary value problem to reformulate the 
problem into a boundary value problem in the unbounded domain in Section 3. 
In Section 4, we present the unique solution of normal reflection for this problem 
when the wedge angle is tt/2. In Section 5, we exhibit the local theory of regular 
reflection-diffraction, introduce a stability criterion to determine state (2) at the 
reflection point on the wedge, and present the von Neumann's detachment and 
sonic conjectures. Then we discuss the role of the potential flow equation in the 
shock reflection-diffraction problem even in the level of the full Euler equations 
in Section 6. Based on the local theory, we reduce the boundary value problem 
into a free boundary problem in the context of potential flow in Section 7. In 
Section 8, we describe a global theory for regular reflection-diffraction for potential 
flow, established in Chen-Feldman [331 1341 135] and Bae-Chen-Feldman [3]. In 
Section 9, we discuss some open problems and new mathematics required for further 
developments, which are also essential for solving multidimensional problems in 
conservation laws and other areas in nonlinear partial differential equations. 

2. Mathematical Formulation I: Initial-Boundary Value Problem 

The full Euler equations for compressible fluids in := R + x IR 2 ,< £ R + := 
(0,oo),x £ E 2 , are of the following form: 



where p is the density, v = (u, v) the fluid velocity, p the pressure, and e the internal 
energy. Two other important thermodynamic variables are the temperature 9 and 
the energy S. The notation a ® b denotes the tensor product of the vectors a and 
b. 

Choosing (p, S) as the independent thermodynamical variables, then the con- 
stitutive relations can be written as (e,p, 6) — (e(p, S),p(p, S),6(p, S)) governed 



(2.1) 




by 



QdS = de + pdr = de ^dp. 

p 2 
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For a polytropic gas, 

R 

(2.2) p=( 7 -l)pe, e = c v d, 7 = 1 + -, 
or equivalcntly, 

(2.3) p = p(p,S)=K P ie s ^, e = e(p,S) = ^-pi- 1 e s ^, 

7- 1 

where R > may be taken to be the universal gas constant divided by the effective 
molecular weight of the particular gas, c v > is the specific heat at constant 
volume, 7 > 1 is the adiabatic exponent, and k > is any constant under scaling. 
When a flow is potential, that is, there is a velocity potential $ such that 

v = V x $, 

then the Euler equations for the flow take the form: 

dtp + div(pV x <i>) = 0, (conservation of mass) 

<„,.!,„ , |2 



where 



9 t $ + -|V X $| 2 + i(p) = B a , (Bernoulli's law) 

p 7 " 1 - 1 

i(p) = when7 > 1, 

7 -l 



especially, i{p) = lnp when 7 = 1, by scaling and Bq is the Bernoulli constant, which 
is usually determined by the boundary conditions if such conditions are prescribed. 
From the second equation in (|2.4p . we have 



(2.5) p(D$) = i-^Bo - (d t $+ i|V x $| 2 )). 

Then system (|2.4[) can be rewritten as the following time-dependent potential flow 
equation of second order: 

(2.6) d t p{D<5>) + V • (p(£>$)V$) = 
with (H3J). 

For a steady solution <i> = y(x), i.e., c? t <t> — 0, we obtain the celebrated steady 
potential flow equation in aerodynamics: 

(2.7) V x • (p(V x $)V x $) = 0. 

In applications in aerodynamics, (|2.4|) or (|2.6|) is used for discontinuous so- 
lutions, and the empirical evidence is that entropy solutions of (|2.4p or (|2.6p are 
fairly close to entropy solutions for (|2.ip . provided the shock strengths are small, 
the curvature of shock fronts is not too large, and the amount of vorticity is small 
in the region of interest. Furthermore, we will show in Section 6 that, for the shock 
reflection-diffraction problem, the Euler equations for potential flow is actually ex- 
act in an important region of the solution (see Theorem 6.1 below). 

Then the problem of shock reflection-diffraction by a wedge can be formulated 
as follows: 

Problem 2.1 (Initial-boundary value problem). Seek a solution of system 
(|2.ip satisfying the initial condition at t = 0: 

(0,0,p ,Po), |x 2 | > xitan^,^! > 0, 

(wi,0,pi,pi), x x < 0; 



(2.8) (w,p,p) 



SHOCK REFLECTION-DIFFRACTION 



5 



and the slip boundary condition along the wedge boundary: 

(2.9) vi/ = 0, 

where v is the exterior unit normal to the wedge boundary, and state (0) and (1) 
satisfy 

/ 01n s (Pi-Po)(Pi ~ Po) Pi (7 + - (7- l)Po 

(2.10) ui = \ , — = - — — — — , pi > po- 



That is, given po,po,pi, and 7 > 1, i/ie other variables U\ and pi are determined 
by (|2.10p . particular, the Mach number Mi = U\/c\ is determined by 

(2.11) M\ = 2( - pl - Po)2 

1 A)((7 + l)Pi - (7- l)Po)' 

where c\ = ^/ JPiJpi is the sonic speed of fluid state (1). 




Figure 1. Initial-boundary value problem 



3. Mathematical Formulation II: Boundary Value Problem 

Notice that the initial-boundary value problem (Problem 2.1) is invariant 
under the self-similar scaling: 

(t, x) — ► (at, ax) for any q^O. 

Therefore, we seek self-similar solutions: 

(v,p,p)(i,x) = (v,p, /})(£, 77), (£,77) = -. 

Then the self-similar solutions are governed by the following system: 
' {p(J)t + (pV) n + 2p = 0, 

(pU 2 +p) s + (pUV)r, + 3pU = 0, 
S (pUV) i + ( P V 2 + p)„ + 3pV = 0, 

. + + + + 2 <^ 2 + t?t' = °' 

where q = VU 2 + V 2 , and (U, V) = (u — £, v — 77) is the pseudo-velocity. 
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The eigenvalues of system (|3.1[) are 



V _UV±cy/ 



Ao = — (repeated), A± = 



2 



u K * y ' 1 u 2 



where c — \J jp/p is the sonic speed. 

When the flow is pseudo-subsonic, i.e., q < c, the eigenvalues A± become com- 
plex and thus the system consists of two transport equations and two nonlinear 
equations of hyperbolic-elliptic mixed type. Therefore, system (|3.1I) is hyperbolic- 
elliptic composite-mixed in general. 

Since the problem is symmetric with respect to the axis rj = 0, it suffices to 
consider the problem in the half-plane 77 > outside the half-wedge: 

A := {£ < 0,17 > 0}U {rj > £taxi0 w , £ > 0}. 

Then the initial-boundary value problem (Problem 2.1) in the (t, x)-coordinatcs 
can be formulated as the following boundary value problem in the self-similar co- 
ordinates (£, rj): 

Problem 3.1 (Boundary value problem in the unbounded domain). 

Seek a solution to system (|3.1[) satisfying the slip boundary condition on the wedge 
boundary and the matching condition on the symmetry line r] = 0: 

{U,V)-u = on 9A = {£<0,?7 = 0}U{£>0,77>£tan(9,4, 

the asymptotic boundary condition as £ 2 + ry 2 — > 00 : 



(U + (,V + r),p,p) 



(0,0,p o ,Po), £ > £o,r) > £tan/ 
(Mi,0,_pi,pi), £<£ o ,?7>0. 




Figure 2 . Boundary value problem in the unbounded domain A 



It is expected that the solutions of Problem 3.1 contain all possible patterns 
of shock reflection-diffraction configurations as observed in physical and numerical 
experiments; cf. EH El Effl EH Il08l fT09j fT29j IT37] and the references cited 
therein. 
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4. Normal Reflection 

The simplest case of the shock reflection-diffraction problem is when the wedge 
angle 9 W is tt/2. In this case, the reflection-diffraction problem simply becomes the 
normal reflection problem, for which the incident shock normally reflects, and the 
reflected shock is also a plane. It can be shown that there exist a unique state 
(p2, P2), Pi > Pi, and a unique location of the reflected shock 



(4.1) 



fc = - 



P1U1 



with 



Ml = 



(p2~Pl)(p2 - Pi) 



such that state (2) 



92- pi y P1P2 

(— £, — 77,^2, P2) is subsonic inside the sonic circle with center 



at the origin and radius C2 = \/jP2/p2, and is supersonic outside the sonic circle 
(see Fig. [3}. That is, in this case, the normal reflection solution is unique. 











reflected 


(2) 




location of 


shock 






incident shock 


(1) 


(-€, -»?>I> a >p 2 ) 







Pl"l 

Pi- Pi 



C — Pl'"l 
>0 Pi- Pa 



Figure 3. Normal reflection solution 



In this case, 
(4.2) Ml = 



2(p 2 - Pi) P2 _ (7 + 1)P2 - (7 - l)pi 



P2 ((7 + l)pi - (7 - 1)P2) ' Pi (7 + I) Pi - (7 - 1)P2 



and ^ = t > 1 is the unique root of 



;i + ^M\y - (2 + ^M x 2 )i + 1 = 0, 



that is 
(4.3) 



p 2 _ 4 + (7 + l)Mi 2 + Mi a/16 + (7 + 1) 2 M2 



Pi 2(2 + ( 7 -l)M 1 2 ) 

In other words, given po,po,pi, and 7 > 1, state (2) = (— £, — rj,p2, P2) is uniquely 
determined through (j2T0| - ([2~TT|) and pTT ]) — <(473 ]) . 



5. Local Theory and von Neumann's Conjectures for Regular 
Reflection-Diffraction Configuration 

For a wedge angle 9 W S (0, 7r/2), different reflection-diffraction patterns may 
occur. Various criteria and conjectures have been proposed for the existence of 
configurations for the patterns (cf. Ben-Dor \6,). One of the most important 
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conjectures made by von Neumann |142L [143 in 1943 is the detachment conjecture, 
which states that the regular reflection-diffraction configuration may exist globally 
whenever the two shock configuration (one is the incident shock and the other the 
reflected shock) exists locally around the point Pq (see Fig. 4). 

The following theorem was rigorously shown in Chang-Chen [18] (also see 
Sheng-Yin [127] . Bleakney-Taub [12], Neumann [142] [143] ). 

Theorem 5.1 (Local theory). There exists 9d = 9d(M s , 7) e (0, ir/2) such that, 
when 9 W £ (9 d ,Tr/2), there are two states (2) = (P 2 a , V£,p%, p%) and V 2 b ,p b 2 , p\) 
such that 

m,V 2 a )\>\(U»,Vi)\ and \(U b ,V b )\ < c b , 

where c\ — yyp^/ p\ is the sonic speed. 

Then the conjecture can be stated as follows: 

The von Neumann's Detachment Conjecture: There exists a global regu- 
lar reflection- diffraction configuration whenever the wedge angle 9 W is in (9d,n/2). 

It is clear that the regular reflection-diffraction configuration is not possible 
without a local two-shock configuration at the reflection point on the wedge, so this 
is the weakest possible criterion. In this case, the local theory indicates that there 
are two possible states for state (2). There had been a long debate to determine 
which one is more physical for the local theory; see Courant-Friedrichs [54] , Ben- 
Dor [5], and the references cited therein. 

Since the reflection-diffraction problem is not a local problem, we take a differ- 
ent point of view that the selection of state (2) should be determined by the global 
features of the problem, more precisely, by the stability of the configuration with 
respect to the wedge angle 9 W , rather than the local features of the problem. 

Stability Criterion to Select the Correct State (2): Since the solution 
is unique when the wedge angle 9 W = 7r/2, it is required that our global regular 
reflection- diffraction configuration should be stable and converge to the unique nor- 
mal reflection solution when 9 W — ► ir/2, provided that such a global configuration 
can be constructed. 

We employ this stability criterion to conclude that our choice for state (2) must 
be (t/f, V 2 ,p 2 , P 2 )- I n general, (U 2 , V 2 ,p 2 , p 2 ) may be supersonic or subsonic. If 
it is supersonic, the propagation speeds are finite and state (2) is completely deter- 
mined by the local information: state (1), state (0), and the location of the point 
Pq. This is, any information from the reflected region, especially the disturbance 
at the corner P3, cannot travel towards the reflection point p). However, if it is 
subsonic, the information can reach Po and interact with it, potentially altering 
the reflection-diffraction type. This argument motivated the second conjecture as 
follows: 

The von Neumann's Sonic Conjecture: There exists a regular reflection- 
diffraction configuration when 9 W € (9 s ,Tt/2) for 9 S > 9d such that \{U 2 , V 2 )\ > c 2 
at P . 

This sonic conjecture is based on the following fact: If state (2) is sonic when 
9 W = 9 S , then \{U 2 ,V 2 )\ > c 2 for any 9 W £ (9 s ,ir/2). This conjecture is stronger 
than the detachment one. In fact, the regime between the angles 9 S and 9d is very 
narrow and is only fractions of a degree apart; see Fig. 5 from Sheng-Yin |127j . 
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Figure 4. Regular reflection-diffraction configuration 



90 



Domain of appearing regular reflection 



ID 



0.1 02 0.3 0.4 0.S 0.6 0.7 0.8 0.9 1 
P 'P 



Figure 5. The von Neumann's sonic criterion vs the detachment 
criterion 9 S > 84 when 7 = f .4. 



6. The Potential Flow Equation 

In this section, we discuss the role of the potential flow equation in the shock 
reflection-diffraction problem for the full Eulcr equations. 

Under the Hodgc-Hclmoltz decomposition (U, V) = Vtp + W with V • W = 0, 
the Euler equations (|34[) become 

(64) V • (pVtp) + 2p + V • (pW) = 0, 

(6.2) V(-\X7tp\ 2 +<p) + -X7p = (Vip + W) ■ X7W + (X7 2 tp + I)W, 

2 p 

(6.3) (Vy? + W) ■ Wuj + (1 + Alp)oj = 0, 

(6.4) (V<y9 + W) ■ VS = 0, 

where u — curl W — curl(t/, V) is the vorticity of the fluid, S = c v \n(pp~~*) is the 
entropy, and the gradient V is with respect to the self-similar variables (£, rf) from 
now on. 

When uj = 0, S = const., and W = on a curve T transverse to the fluid 
direction, we first conclude from (|6.3[) that, in the domain fij determined by the 
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fluid trajectories: 
past r, 

uj = 0, i.e. curl IF = 0. 
This implies that W — const, since V • W = 0. Then we conclude that 

W = infix, 

since TF|r = 0, which yields that the right-hand side of equation (|6.2|) vanishes. 
Furthermore, from (|6.4p . 

S = const. in fix, 

which implies that 

p = const, p 1 . 

By scaling, we finally conclude that the solutions of system (|6.ip - (|6.4[) in the domain 
fix is determined by the following system for self-similar solutions: 

f V ■ (i>V V ) + 2p = 0, 

(6 - 5) 

or the potential flow equation for self-similar solutions: 

(6.6) V-(p{V<p,p)V<p)+2p{V<p,p) = 0, 
with 

(6.7) P (\Vp\ 2 ,p) = {pr 1 - { 1 -l){i P + l -\V^))^. 
Then we have 

(6.8) c 2 = c 2 (|V^| 2 , p, p^ 1 ) = pr 1 (7 - 1)(||V^| 2 + <p). 
For our problem (see Fig. 4), we note that, for state (2), 

(6.9) lo = 0, W = 0, S = S 2 - 

Then, if our solution (U, V, p, p) is C ' 1 and the gradient of the tangential component 



of the velocity is continuous across the sonic arc r so „i C , we still have f|(3. 9|) along 
r 'sonic on the side of fi. Thus, we have 

Theorem 6.1. Let (U,V,p,p) be a solution of our Problem 3.1 such that 
(U,V,p,p) is C 0,1 in the open region P0P1P2P3 and the gradient of the tangen- 
tial component of (U, V) is continuous across the sonic arc T son i C . Let fix be the 
subregion of fi formed by the fluid trajectories past the sonic arc T 'sonic, then, in 
fix, the potential flow equation (|6.6p with (|6.7[) coincides with the full Euler equa- 
tions (|6.ip - (|6.4|) , that is, equation (|6 . 6[) with (|6.7|) is exact in the domain fix for 
Problem 3.1. 

Remark 6.1. The regions such as fix also exist in various Mach reflection- 
diffraction configurations. Theorem 16.11 applies to such regions whenever the solu- 
tion (£/, V, p, p) is C 0,1 and the gradient of the tangential component of (U,V) is 
continuous. In fact, Theorem 8.3 indicates that, for the solutions tp of (|6.6|) with 
(|6.7p . the C ' regularity of ip and the continuity of the tangential component of 
the velocity field (U, V) — Vp are optimal across the sonic arc T SO nic- 
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Remark 6.2. The importance of the potential flow equation (12. 6|) with (J275J) 
in the time-dependent Euler flows was also observed by Hadamard [83j through a 
different argument. 

Furthermore, when the wedge angle 9 W is close to 7r/2, it is expected that 
the curvature of the reflected shock is small so that, in the other part of f2, 
the vorticity to is small and the entropy is close to the constant. Then, in the 
reflection-diffraction domain Q = 1]] U U2, the potential flow equation (|6.6|) with 
()6.7p dominates, provided that the exact state along the reflected shock is given. 




Potential Flow 



Figure 6. The potential flow equation dominates the domain Q 

Equation (|6.6[) with (|6.7p is a nonlinear equation of mixed elliptic-hyperbolic 
type. It is elliptic if and only if 

(6.10) \V<p\<c{\V<p\ 2 ,<p,f%- 1 ), 
which is equivalent to 



(6.11) \V<p\ < c,(^po,7) := 




-( 7 -l)p). 



The study of partial differential equations of mixed hyperbolic-elliptic type can 
date back 1940s (cf. pH HI IT461 [149] ). Linear models of partial differ- 
ential equations of mixed hyperbolic-elliptic type include the Lavrentyev-Betsadze 
equation: 

d xx u + sign(x)dyyii = 0, 

the Tricomi equation: 

u xx + xuyy = (hyperbolic degeneracy at x = 0), 

and the Keldysh equation: 

xu xx + u yy = (parabolic degeneracy at x = 0) . 

Nonlinear models of mixed-type equations for (|6.6p with (|6.7p include the transonic 
small disturbance equation: 




or, for v — u — x, 

, ^ 3 1 n 

(V V X ) X + Vyy + ~V X + ~ = 0, 
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which has been studied in [HI US [87), [83 \89\ I116j and the references cited 
therein. Also see |16L I130|, I157|, I158j for the models for self-similar solutions from 
the pressure gradient system and nonlinear wave equations. 

7. Mathematical Formulation III: Free Boundary Problem for 

Potential Flow 



For the potential equation (|6.6|) with (|6.7|) . shocks are discontinuities in the 
pseudo-velocity Vip. That is, if D + and D~ := D \ D + are two nonempty open 
subsets of D C ffi 2 and S := dD + n D is a C 1 -curve where Dip has a jump, then 
if E Wj£(D) n C 1 {D ± US)D C 2 (D ± ) is a global weak solution of (EU) with (15771) 
in Z) if and only if 99 is in W; '^°(£>) and satisfies equation ()6.6p in D^ 1 and the 
Rankine-Hugoniot condition on 5: 

(7.1) [p(.\V<p\ 2 ,<p)V<P-"] s =0, 

where the bracket [■] denotes the difference of the values of the quantity along the 
two sides of S. 

Then the plane incident shock solution in the (t, x)-coordinates with states 
(p, V x $) = (po, 0, 0) and (pi, u\, 0) corresponds to a continuous weak solution p of 
(16.5|) in the self-similar coordinates (£, 77) with the following form: 

(7.2) Mtv)=-\(e+V 2 ) for 

(7.3) <Pi(d,v)=-\(e+V 2 )+Mt-to) for £<£o, 
respectively, where 



(7.4) "^ f'Vlfl'f^ - 

V (7-l)(Pi+A)) 



(7.5) fr=pw ;r\w/° 2 ; = ^>o 



2 (Pi 1 - Po ]) = 
(7-l)(p?-Po) P1-P0 
are the velocity of state (1) and the location of the incident shock, uniquely de- 
termined by (po,pi,j) through (|7.1[) . Then Po = (£o> £0 tan6*„,) in Fig. 2, and 
Problem 3.1 in the context of the potential flow equation can be formulated as: 

Problem 7.1 (Boundary value problem) (see Fig. 2). Seek a solution ip of 
equation (|6 - 6[) with (|6.7p in the self-similar domain A with the boundary condition 
on OA: 

(7.6) V<p-v\ 8 A = 0, 

and the asymptotic boundary condition at infinity: 

- j<A) for i > £o,i] > ^tan6» w , 2 2 

[Vi /or 4 < Co,?? > 0, 
where (|7.7[) ZioZrfs m £/ie sense that lim — ^||c(A\B K (o)) = 0- 

For our problem, since <p>\ does not satisfy the slip boundary condition (|7.6p . 
the solution must differ from epi in {£ < Co} HA, thus a shock diffraction-diffraction 
by the wedge vertex occurs. In Chen-Feldman [33|, 134) , we first followed the von 
Neumann criterion and the stability criterion introduced in Section 5 to establish 
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a local existence theory of regular shock reflection near the reflection point Pq in 
the level of potential flow, when the wedge angle is large and close to ir/2. In this 
case, the vertical line is the incident shock S = {£ = £0} that hits the wedge at the 
point Pq = (£0, £0 tan 9 W ), and state (0) and state (1) ahead of and behind S are 
given by tpo and tp\ defined in (|7.2|) and (|7.3p . respectively. The solutions tp and 
(pi differ only in the domain P0P1P2P3 because of shock diffraction by the wedge 
vertex, where the curve P0P1P2 is the reflected shock with the straight segment 
PqPi- State (2) behind P3P1 can be computed explicitly with the form: 

(7.8) tp 2 (£, V) = ~\(? + V 2 ) + - Co) + (V - £0 tand„)«a tan6 w , 

which satisfies Vtp ■ v = on <9A n {£, > 0}; the constant velocity 112 and the 
angle 9 S between PqP\ and the £-axis are determined by (6 W , po , pi , 7) from the 
two algebraic equations expressing (|7.I|) and continuous matching of state (1) and 
state (2) across PqP\, whose existence is exactly guaranteed by the condition on 
(0w: Pot Pi: 7) under which regular shock reflection-diffraction is expected to occur 
as in Theorem 5.1. Moreover, (p% is the unique solution in the domain PqP±P4, as 
argued in [181 1123] . Denote 

P1P4 := r ' son i C =dttn dB C2 (u-2, u 2 tan 6 W ) 

the sonic arc of state (2) with center (1*2, 1(2 tan W ) and radius C2. Also we introduce 
the following notation for the other parts of dCl: 

Tshock ■= P1P2; r wedge : = dnndAn{ v >o} = p 3 p 4 ; r symm ■.= { v = o}ndn. 

Then Problem 7.1 can be formulated as: 

Problem 7.2. Seek a solution tp in Q to equation (|6.6|) with (|6.7p subject to 
the boundary condition (|7.6|) on dfl n dA, the Rankine-Hugoniot conditions on the 

shock r shock ■' 

(7-9) [<p]r, hock - 0, 

(7-10) WV^^pojVr^^O, 
and the Dirichlet boundary condition on the sonic arc T son i c : 
(7.11) (¥>-¥>2)|r so „ io =0. 



It should be noted that, in order that the solution ip in the domain fl is a part 
of the global solution to Problem 7.1, that is, ip satisfies the equation in the sense 
of distributions in A, especially across the sonic arc T son i Cl it is requaired that 

V(ip - ip 2 ) ■ v\r sonlc = 0. 

That is, we have to match our solution with state (2), which is the necessary 
condition for our solution in the domain f2 to be a part of the global solution. To 
achieve this, we have to show that our solution is at least C 1 with V(y — P2) = 

aCrOSS r sonic- 

Then the problem can be reformulated as the following free boundary problem: 

Problem 7.3 (Free boundary problem). Seek a solution tp and a free 
boundary T s i loc k — {£ = /(??)} such that 

(i) / e C 1 '" and 
(7.12) n+ = {£ >f(r,)} HD={p<p 1 }n D; 
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(ii) tp satisfies the free boundary condition (|7.10|) along T 'shock ; 

(iii) ip g C 1 ' Q (r2 + )nC 2 (f2_|_) solves (|6.5p mf2 + , is subsonic infl + , and satisfies 

(7-13) {<P-<PiM<P-<P2)-v)\v.^=0, 
(7-14) • Hrw ge ur symm = 0. 

The boundary condition on T symrn implies that /'(0) — and thus ensures the 
orthogonality of the free boundary with the £-axis. Formulation (|T. 1 2[) implies that 
the free boundary is determined by the level set ip = ipi, which is a convenient 
formulation to apply useful free boundary techniques. The free boundary condition 
(|7.10|) along T s hock is the conormal boundary condition on T shock- Condition (|7.13j) 
ensures that the solution of the free boundary problem in £1 is a part of the global 
solution as pointed out earlier. Condition (IT. 14[) is the slip boundary condition. 

Problem 7.3 involves two types of transonic flow: one is a continuous transi- 
tion through the sonic arc r 'sonic as a fixed boundary from the pseudo-supersonic 
region (2) to the pseudo-subsonic region f); the other is a jump transition through 
the transonic shock as a free boundary from the supersonic region (1) to the sub- 
sonic region Q. 

8. Global Theory for Regular Reflection-Diffraction for Potential Flow 

In this section, we describe a global theory for regular reflection-diffraction 
established in Chen-Feldman [33 |, I34| , [35] and Bae-Chen-Feldman [3]. 

8.1. Existence and stability of regular reflection-diffraction configu- 
rations. In Chen-Feldman [33|, I34j . we have developed a rigorous mathematical 
approach to solve Problem 7.3 and established a global theory for solutions of 
regular reflection-diffraction, which converge to the unique solution of the normal 
shock reflection when 9 W tends to ir/2. 

Introduce the polar coordinates (r, 9) with respect to the center (u2, U2 ta,rx6 w ) 
of the sonic arc Y sonic of state (2), that is, 

(8.1) £ — it2 = rcos#, rj — U2ta,Ti9 w = rsm.9. 
Then, for e € (0,C2), we denote by 

Q £ :=nn{{r,9) : < c 2 - r < e} 

the e-neighborhood of the sonic arc P1P4 within fi; see Fig. 4. In fl £ , we introduce 
the coordinates: 

(8.2) x = c 2 -r, y = 6- W . 

Then fl E C {0 < x < e, y > 0} and P1P4 C {x = y > 0}. 

Theorem 8.1 (Chen-Feldman [33l [34]). There exist C = 9 c (p ,Pi,l) € 
(0,7r/2) and a = a(po, pi,j) € (0, 1/2) such that, when 9 W € (6 c ,tt/2), there exists 
a global self-similar solution 

$(i,x) =t<p(-) + ^- for -e A, <>0 
y ; ^ K t' 2t t 

with p(t,x.) — (pq 1 ~ — ^l^x^l 2 )^ 1 " o/Problem 7.1 (equivalently, Problem 
7.2) for shock reflection- diffraction by the wedge. The solution tp satisfies that, for 
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(8.3) 



(&»?) = x/t, 

<p G C° :1 (A), 

{'Po / or £ > £o an d V > 6 tan , 
/ or £ < Co arl ^ above the reflected shock P0P1P2, 
<p 2 in PqPiP 4 - 

Moreover, 

(i) equation (16. 5|) is elliptic in fi; 

(ii) </?2 < y < Vi i n 

(iii) </ie reflected shock PqP\P2 is C 2 at P\ and C°° except P\; 



(iv) there exists £q 6 (0, C2/2) suc/i i/iai 6 C ' (fl Eo ) H C (Q £g \ r s0 mc)/ in 
particular, in the coordinates 



30; 



(8.4) 1^-^11^:= J2 SU P (^ + ^ 2 ^>-^)(z,y)|) < 

(v) there exists 5q > so that, in the coordinates (|8.2[) , 

2 

(8.5) \d x (ip - ip 2 )(x,y)\ < — — -x in n E0 ; 

7+1 

(vi) there exist u) > and a function y = f{x) such that, in the coordinates 

(El, 



(8.6) 



O E0 = {(x,y) : x € (0, e„), < j/ < />)}, 
r s / iOC /c n dn eo = {(x,y) : x e (0, e ), y = /(x)}, 



df 

.7) ||/||ci.i([o, eo]) < °°> — > w > /or < X < £ - 



Furthermore, the solution (p is stable with respect to the wedge angle 9 W in W, 1 



AA 

' loc 

and converges in Wj ' c to the unique solution of the normal reflection as 9 W — » 7r/2. 

The existence of a solution ip of Problem 7.1 (equivalently, Problem 7.2), 
satisfying (|8 . 3[) and property jlv|, follows from [34i Main Theorem]. Property (0) 
follows from Lemma 5.2 and Proposition 7.1 in 34 . Property (JTTJ) follows from 
Proposition 7.1 and Section 9 in [34] which assert that where the set 

/C is defined by (5.15) in [34] . Property (jvj follows from Propositions 8.1-8.2 and 
Section 9 in [34]. Property flvl]) follows from (5.7) and (5.25)-(5.27) in [34] and the 
fact that cp — (fi2 £ IC. 

We remark that estimate ([8.4]) above confirms that our solutions satisfy the 
assumptions of Theorem 6.1 for the velocity field (U, V) = "Vip. 

One of the main difficulties for the global existence is that the ellipticity con- 
dition ([6.10]) for (]6.6|) with (|6.7p is hard to control, in comparison to our work on 
steady flow [29], [30l ETQ [32] . The second difficulty is that the ellipticity degener- 
ates along the sonic arc T son i C . The third difficulty is that, on T scm i C , the solution 
in has to be matched with ip 2 at least in C , i.e., the two conditions on the 
fixed boundary r so „i C : the Dirichlet and conormal conditions, which are generi- 
cally overdetermined for an elliptic equation since the conditions on the other parts 
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of boundary have been prescribed. Thus, one needs to prove that, if ip satisfies 
(|6.6p in Q, the Dirichlct continuity condition on the sonic arc, and the appropriate 
conditions on the other parts of dfl derived from Problem 7.3, then the normal de- 
rivative • v automatically matches with V</?2 ■ v along T SO nic- Indeed, equation 
(|6.6p . written in terms of the function tp — ip — ipi i n the (x, y)-coordinates defined 
near I" sonic such that T son ic becomes a segment on {x = 0}, has the form: 



plus the small terms that are controlled by n/2—d w in appropriate norms. Equation 
(18. 8[) is elliptic if tp x < 2x/(7+ 1). Hence, it is required to obtain the C 1 ' 1 estimates 
near T SO nic to ensure \?p x \ < 2x/(7 + 1) which in turn implies both the ellipticity 
of the equation in Q and the match of normal derivatives V<p • v = V</?2 ■ v along 
r sonic- Taking into account the small terms to be added to equation (|8.8|) . one needs 
to make the stronger estimate \ip x \ < 4a;/ (3(7 + 1)) and assume that tt/2 — 9 W is 
suitably small to control these additional terms. Another issue is the non- variational 
structure and nonlinearity of this problem which makes it hard to apply directly 
the approaches of Caffarelli [13j and Alt-Caffarelli- Friedman [H [2]. Moreover, the 
elliptic degeneracy and geometry of the problem makes it difficult to apply the 
hodograph transform approach in Chen-Feldman [3 1 1 and Kinderlehrer-Nirenberg 
|94) to fix the free boundary. 

For these reasons, one of the new ingredients in our approach is to develop 
further the iteration scheme in [291 131] to a partially modified equation. We 
modified equation (|6 . 6[) in Q by a proper Shiffmanization (i.e. a cutoff) that depends 
on the distance to the sonic arc, so that the original and modified equations coincide 
when <p satisfies \ip x \ < 4a;/ (3(7+ 1)), and the modified equation Afcp = is elliptic 
in Q with elliptic degeneracy on P1P4. Then we solved a free boundary problem 
for this modified equation: The free boundary is the curve T shock, and the free 
boundary conditions on T s hock are ip = ip\ and the Rankine-Hugoniot condition 



On each step, an iteration free boundary curve r so „i C is given, and a solution of 
the modified equation Afip = is constructed in fl with the boundary condition (|7.ip 
on r shock, the Dirichlet condition tp = ip2 on the degenerate arc T SO nic, and Vtp-v — 
on P2P3 and T we d ge - Then we proved that ip is in fact C ' up to the boundary 
r sonic, especially \V(<p — <y3 2 )| < Cx, by using the nonlinear structure of elliptic 
degeneracy near T son ic which is modeled by equation (|8.8p and a scaling technique 
similar to Daskalopoulos-Hamilton 56.' and Lin- Wang [105j . Furthermore, we 
modified the iteration free boundary curve T shock by using the Dirichlet condition 
<p = (pi on r shock ■ A fixed point tp of this iteration procedure is a solution of the 
free boundary problem for the modified equation. Moreover, we proved the precise 
gradient estimate: \ip x \ < 4a;/ (3(7 + 1)) for which implies that <p satisfies the 
original equation (|6.5p . 

This global theory for large-angle wedges has been extended in Chen-Feldman 
[35j to the sonic angle 6 S < 9 C , for which state (2) is sonic, such that, as long as 
6 W £ (9 s , tt/2], the global regular reflection-diffraction configuration exists. 



(8 




in x > and near x = 0, 



Theorem 8.2 (von Neumann's Sonic Conjecture (Chen-Feldman [35])). The 
global existence result in Theorem 8.1 can be extended up to the sonic wedge-angle 
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9 a for any 7 > 1 and u\ < c\. Moreover, the solutions satisfy the properties (i)-(vi) 
in Theorem 8.1. 

The condition u\ < c\ depends explicitly only on the parameters 7 > 1 and 
pi > pa > 0. For the case u\ > c\, we have been making substantial progress as 
well, and the final detailed results can be found in Chen-Fcldman [35 . 

8.2. Optimal regularity. By Theorem 8.1([jv]), the solution ip constructed 
there is at least C ' near the sonic arc T son i C . The next question is to analyze the 
behavior of solutions <p{£,,n) to regular reflection-diffraction, especially the optimal 
regularity of the solutions. 

We first define the class of regular reflection-diffraction solutions. 

Definition 8.3. Lei 7 > 1, pi > p > 0, andO w e (0,7r/2) be constants, let ui 
and £0 be defined by (|7.4p and (|7.5|) . Let the incident shock S = {£ = £0} hits the 
wedge at the point Pq = (£cb£o ta,n9 w ), and let state (0) and state (1) ahead of and 
behind r s ^ oc ^. be given by (|7.2|) and (|7.3|) . respectively. The function ip £ C 0,1 (A) 
is a regular reflection- diffraction solution if ip is a solution to Problem 7.1 such 
that 

(a) there exists state (2) of form (|7.8p with U2 > 0, satisfying the entropy 
condition P2 > pi and the Rankine-Hugoniot condition (|7.1|) along the line Si := 
{pi = 1P2} which contains the points Pq and Pi, such that Pi G A is on the sonic 
circle of state (2), and state (2) is supersonic along PqPi; 

(b) there exists an open, connected domain Q := P1P2P3P4 C A such that (|8.3p 
holds and equation (|6.5[) is elliptic in J7; 

(c) p > P2 on the part P1P2 = T shock of the reflected shock. 

Remark 8.1. The global solution constructed in |33L 1341 135) is a regular 
reflection-diffraction solution, which is a part of the assertions in Theorems 8.1- 
8.2. 

Remark 8.2. If state (2) exists and is supersonic, then the line Si = {tpi = P2] 
necessarily intersects the sonic circle of state (2); see the argument in [331 134] 
starting from (3.5) there. Thus, the only assumption regarding the point P\ is that 
Si intersects the sonic circle within A. 

Remark 8.3. We note that, in the case 9 W — the regular reflection becomes 
the normal reflection, in which U2 = and the solution is smooth across the sonic 
line of state (2); see |34l Section 3.1]. Condition 9 W £ (0, in Definition 18.31 rules 
out this case. Moreover, for 9 W € (0, f), the property 112 > in part (a) of Definition 
I8.3l is always true for state (2) of form (|7.8|) . satisfying the entropy condition p 2 > p\ 
and the Rankine-Hugoniot condition (7.1) along the line S\ := {ipi = ^2} which 
contains the point Pq. These are readily derived from the calculations in [341 
Section 3.2]. 

Remark 8.4. There may exist a global regular reflection-diffraction configura- 
tion when state (2) is subsonic which is a very narrow regime |54L \142\ [143] . Such 
a case does not involve the difficulty of elliptic degeneracy, which we are facing for 
the configurations in the class of solutions in the sense of Definition 8.3. 

Remark 8.5. Since ip = ip% on T s h oc k by (|8.3|) . condition (c) in Definition 18.31 
is equivalent to 

^shock C {ip 2 < <ySl}, 
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that is, r shock is below Si. 

In Bae-Chen-Feldman [3], we have developed a mathematical approach to es- 
tablish the regularity of solutions of the regular reflection-diffraction problem in the 
sense of Definition 8.3. 

First, we have shown that any regular reflection-diffraction solutions cannot be 
C 2 across the sonic arc r sonic := P1P4. 

Theorem 8.4 (Bae-Chcn-Fcldman [3]). There does not exist a global regular 
reflection- diffraction solution in the sense of Definition 8.3 such that tp is C 2 across 
the sonic arc T SO nic- 

Now we study the one-sided regularity up to r so „i C from the elliptic side, i.e., 
from Q. For simplicity of presentation, we now use a localized version of For a 
given neighborhood N"(T SO nic) of T son ic and e > 0, define 

n E := fi njV(r sonic ) n {x < e}. 

Since Af(T son i c ) will be fixed in the following theorem, we do not specify the de- 
pendence of fl £ On J\f(T gonic). 

Theorem 8.5 (Bae-Chen-Feldman [3]). Let ip be a regular reflection- diffraction 
solution in the sense of Definition 18.31 and satisfy the following properties: There 
exists a neighborhood Af (T SO nic) of T sonic such that 

(a) ip is C ' across the sonic arc T son i C : ip G C ' (P0P1P2P3 H Al '(T sonic)) ; 

(b) there exists Sq > so that, in the coordinates (|8.2|) . 

2 S 

(8.9) \d x (<p - if 2 )(x,y)\ < — — -x in nnAf(T sonic); 

7 + 1 

(c) there exists eg > 0, lu > 0, and a function y = f(x) such that, in the 
coordinates (|8.2p . 



(8.10) 



n eo = {(x,y) : x € (0, e ), < y < f(x)}, 
Fshock n dVL eo = {(x,y) : x G (0, e ), V = /(^)}, 



and 



1 r 

(8.11) ||/|| c i,i ([0 , £o]) < 00, — > u > for < x < e . 

Then we have 

(i) ip is C 2 a up to Tsonic away from P\ for any a G (0, 1). That is, for any 
a G (0, 1) and any given (£o)Vo) G T son ic \ {Pi}, there exists K < 00 de- 
pending only on p , pi, 7, e , a, W^Wc 1 - 1 ^), andd = dist((£ ,rj ), T s hock) 
so that 

(ii) For any (€o,Vo) G r sonic \ {PJ, 

1 



lim (D rr (p — D rr p>2) 



K.>7)-'«0.')0) 7 + 1 ' 
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(iii) D 2 ip has a jump across T sonic : For any (cjo,?7o) € T sonic \ {Pi}, 

lim D rr ip — lim D rr ip = . 

((,v)->(So,io) (t,v)^(iO,Vo) 7 + 1' 

(£,»7)sn (€,<7)eA\f! 

lim (D r9 , D eg )(p = lim (D r6 ,D ee )ip = 0; 

(£,»j)en (5,i7)eA\s-2 

(iv) TTie Zimitlimte^j^P! _D 2 (/? does not exist. 

We remark that the solutions established in [331 1341 [35] satisfy the assumptions 
of Theorem 8.5. In particular, we proved that the C 1: ^regularity is optimal for the 
solution across the open part P1P4 of the sonic arc (the degenerate elliptic curve) 
and at the point P\ where the sonic circle meets the reflected shock (as a free 
boundary). 

To achieve the optimal regularity, one of the main difficulties is that the sonic 
arc r 'sonic is the transonic boundary separating the elliptic region from the hy- 
perbolic region, near where the solution is governed by the nonlinear degenerate 
elliptic equation (|8.8p for tp = (p — (f2. We carefully analyzed the features of equa- 
tion (|8.8[) and established the C 2,a regularity of solutions in the elliptic region up 
to the open sonic arc P1P4. As a corollary, we showed that the C 1 ' 1 -regularity is 
actually optimal across the transonic boundary P1P2 from the elliptic to hyperbolic 
region. Since the reflected shock P\Pi is regarded as a free boundary connecting 
the hyperbolic region (1) with the elliptic region fl for the nonlinear second-order 
equation of mixed type, another difficulty for the optimal regularity of the solution 
is that the point Pi is exactly the point where the degenerate elliptic arc P1P4 meets 
a transonic free boundary for the nonlinear partial differential equation of second 
order. As far as we know, this is the first optimal regularity result for solutions 
to a free boundary problem of nonlinear degenerate elliptic equations at the point 
where an elliptic degenerate curve meets the free boundary. To achieve this, we 
carefully constructed two sequences of points on where the corresponding sequences 
of values of ip xx have different limits at Pi; this has been done by employing the 
one-sided C 2,a regularity of the solution up to the open arc P1P4 and by studying 
detailed features of the free boundary conditions on the free boundary P1P2, i.e., 
the Rankine-Hugoniot conditions. 

We remark that some efforts were also made mathematically for the reflection- 
diffraction problem via simplified models. One of these models, the unsteady tran- 
sonic small-disturbance (UTSD) equation, was derived and used in Keller-Blank 
[93) . Hunter-Keller [90] , Hunter [87] . and Morawetz [116] for asymptotic analy- 
sis of shock reflection-diffraction. Also see Zheng |158j for the pressure gradient 
equation and Canic-Keyfitz-Kim [15] for the UTSD equation and the nonlinear 
wave system. Furthermore, in order to deal with the reflection-diffraction prob- 
lem, some asymptotic methods have been also developed. Lighthill [104] studied 
shock reflection-diffraction under the assumption that the wedge angle is either 
very small or close to tt/2. Keller-Blank [93] , Hunter-Keller [90] . Harabetian [84] . 
and Gamba-Rosales-Tabak |67] considered the problem under the assumption that 
the shock is so weak that its motion can be approximated by an acoustic wave. 
For a weak incident shock and a wedge with small angle for potential flow, by 
taking the jump of the incident shock as a small parameter, the nature of the 
shock reflection-diffraction pattern was explored in Morawetz [116 by a number 
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of different scalings, a study of mixed equations, and matching the asymptotics 
for the different scalings. Also see Chen [39 for a linear approximation of shock 
reflection-diffraction when the wedge angle is close to 7r/2 and Serre [123j for an 
apriori analysis of solutions of shock reflection-diffraction and related discussions 
in the context of the Euler equations for isentropic and adiabatic fluids. 

Another related recent effort has been on various important physical problems 
in steady potential flow, as well as steady fully Euler flow, for which great progress 
has been made. The problems for global subsonic flow past an obstacle and for 
local supersonic flow past an obstacle with sharp head are classical, due to the 
works of Shiffman 128], Bers [9], Finn-Gilbarg |64| , Dong [59] . and others for the 
pure elliptic case and to the works of Gu [79] , Shaeffer [119j . Li [102] , and others 
for the pure hyperbolic case. Recent progress has been made on transonic flow past 
nozzles (e.g. HH [26l US |30l [32l H3 SH SS [148], fT50j H51] ). transonic flow 
past a wedge or conical body (e.g. [28L 1471 163]). and transonic flow past a smooth 
obstacle (e.g. [36], [66l fTT5l HT7] V Also see gTJ HIT] for the existence of global 
subsonic-sonic flow, [38l H3 ED IM HI [48], HUB"! IT55] for global supersonic flow 
past an obstacle with sharp head, and the references cited therein. 

For some of other recent related developments, we refer the reader to Chen 
|45L 146] for the local stability of Mach configuration, Elling-Liu 162] for physical- 
ity of weak Prandtl-Meyer reflection for supersonic potential flow around a ramp, 
Serre [125j for multidimensional shock interaction for a Chaplygin gas, Canic- 
Keyfitz-Kim |15L 116] for semi-global solutions for the shock reflection problem, 
Glimm-Ji-Li-Li-Zhang-Zhang-Zheng [75j for the formation of a transonic shock 
in a rarefaction Riemann problem for polytropic gases, Zheng [156] 1157] 1158] 
for various solutions to some two-dimensional Riemann problems, Gues-Metivier- 
Williams-Zumbrun [811 182] and Benzoni-Gavage and Serre [7] for revisits of the 
local stability of multidimensional shocks and phase boundaries, among many oth- 
ers. 

9. Shock Reflection-Diffraction vs New Mathematics 

As we have seen from the previous discussion, the shock reflection-diffraction 
problem involves several core challenging difficulties that we have to face in the 
study of nonlinear partial differential equations. These nonlinear difficulties in- 
clude free boundary problems, oblique derivative problems for nonsmooth domains, 
degenerate elliptic equations, degenerate hyperbolic equations, transport equations 
with rough coefficients, mixed and/or composite equations of hyperbolic- elliptic type, 
behavior of solutions when a free boundary meets an elliptic degenerate curve, and 
compressible vortex sheets. More efficient numerical methods are also required for 
further understanding of shock reflection-diffraction phenomena. 

Furthermore, the wave patterns of shock reflection-diffraction configurations are 
the core patterns and configurations for the global solutions of the two-dimensional 
Riemann problem; these solutions are building blocks and local structure of gen- 
eral entropy solutions and determine global attractors and asymptotic states of 
entropy solutions, as time goes infinity, for two-dimensional systems of hyperbolic 
conservation laws. 

Therefore, a successful solution to the shock reflection-diffraction problem not 
only provides our understanding of shock reflection-diffraction phenomena and be- 
havior of entropy solutions to multidimensional conservation laws, but also provides 
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important new ideas, insights, techniques, and approaches for our developments of 
more efficient analytical techniques and methods to overcome the core challenging 
difficulties in multidimensional problems in conservation laws and other areas in 
nonlinear partial differential equations. The shock reflection-diffraction problem is 
also an excellent test problem to examine our capacity and ability to solve rigor- 
ously various challenging problems for nonlinear partial differential equations and 
related applications. 



References 

[I] H. W. Alt and L. A. CafTarelli, Existence and regularity for a minimum problem with free 
boundary, J. Reine Angew. Math. 325 (1981), 105-144. 

[2] H. W. Alt, L. A. CafTarelli, and A. Friedman, A free-boundary problem for quasilinear elliptic 

equations, Ann. Scuola Norm. Sup. Pisa CI. Sci. (4), 11 (1984), 1-44. 
[3] M.-J. Bae, G.-Q. Chen, and M. Feldman, Regularity of solutions to regular shock reflection for 

potential flow, Inventiones Mathematicae, 2009 (to appear); DOI 10.1007/s00222-008-0156-4 

(online); Preprint larXiv:0804.2500l April 2008. 
[4] M.-J. Bae and M. Feldman, Transonic shocks of multidimensional compressible flow through 

divergent nozzles with arbitrary cross-sections, In preparation, 2008. 
[5] G. Ben-Dor,, Shock Wave Reflection Phenomena, Springer- Verlag: New York, 1991. 
[6] T. V. Bazhcnova, L. G. Gvozdeva, and M. A. Nettleton, Unsteady interaction of shock waves, 

Prog. Aero. Sci. 21 (1984), 249-331. 
[7] S. Benzoni-Gavage and D. Serre, Multidimensional Hyperbolic Partial Differential Equations. 

First-Order Systems and Applications, The Clarendon Press and Oxford University Press: 

Oxford, 2007. 

[8] M. Berger and P. Colella, Local adaptive mesh refinement for shock hydrodynamics, J. Comp. 
Phys. 82 (1989), 64-84. 

[9] L. Bers, Existence and uniqueness of subsonic flows past a given profile, Comm. Pure Appl. 

Math. 7 (1954), 441-504. 
[10] L. Bers, Mathematical Aspects of Subsonic and Transonic Gas Dynamics, John Wiley & 

Sons, Inc.: New York; Chapman & Hall, Ltd.: London 1958. 

[II] A. V. Betsadze, Equations of the Mixed Type, Macmillan Company: New York, 1964. 

[12] W. Bleakney and A. H. Taub, Interaction of shock waves, Rev. Modern Phys. 21 (1949), 
584-605. 

[13] L. A. CafTarelli, A Harnack inequality approach to the regularity of free boundaries, I. Lip- 
schitz free boundaries are C 1,Q! , Rev. Mat. Iberoamericana, 3 (1987), 139-162; II. Flat free 
boundaries are Lipschitz, Comm. Pure Appl. Math. 42 (1989), 55-78; Existence theory, 
compactness, and dependence on X, Ann. Scuola Norm. Sup. Pisa CI. Sci. (4), 15 (1989), 
583-602. 

[14] L. A. CafTrelli, D. Jerison, and C. Kenig, Some new monotonicity theorems with applications 

to free boundary problems, Ann. Math. (2) 155 (2002), 369-404. 
[15] S. Canic, B. L. Keyfitz, and E. H. Kim, A free boundary problems for a quasilinear degenerate 

elliptic equation: regular reflection of weak shocks, Comm. Pure Appl. Math. 55 (2002), 71— 

92. 

[16] S. Canic, B. L. Keyfitz, and E. H. Kim, Free boundary problems for nonlinear wave systems: 
Mach stems for interacting shocks, SIAM J. Math. Anal. 37 (2006), 1947-1977. 

[17] Canic, S., Keyfitz, B. L., and Lieberman, C, A proof of existence of perturbed steady tran- 
sonic shocks via a free boundary problem, Comm. Pure Appl. Math. 53 (2000), 484-511. 

[18] T. Chang, and G.-Q. Chen, Diffraction of planar shock along the compressive corner, Acta 
Math. Scientia, 6 (1986), 241-257. 

[19] T. Chang, G.-Q. Chen, and S. Yang, 2-D Riemann problem in gas dynamics and forma- 
tion of spiral, In: Nonlinear Problems in Engineering and Science-Numerical and Analytical 
Approach (Beijing, 1991), 167-179, Science Press: Beijing, 1992. 

[20] T. Chang, G.-Q. Chen, and S. Yang, On the Riemann Problem for two-dimensional Euler 
equations I: Interaction of shocks and rarefaction waves; II: Interaction of contact disconti- 
nuities, Discrete Contin. Dynam. Systems, 1 (1995), 555-584; 6 (2000), 419-430. 



22 GUI-QIANG CHEN AND MIKHAIL FELDMAN 

[21] T. Chang and L. Hsiao, The Riemann Problem and Interaction of Waves in Gas Dynamics, 
Longman Scientific & Technical, Harlow; and John Wiley & Sons, Inc.: New York, 1989. 

[22] G.-Q. Chen, Euler-Poisson-Darboux equations and hyperbolic conservation laws, In: Non- 
linear Evolutionary Partial Differential Equations, pp. 11-25, AMS/IP Studies in Advanced 
Mathematics, Vol. 3, International Press: Cambridge, Massachusetts, 1997. 

[23] G.-Q. Chen, Euler Equations and Related Hyperbolic Conservation Laws, Chapter 1, Hand- 
book of Differential Equations, Evolutionary Equations, Vol. 2, Eds. C. M. Dafcrmos and E. 
Feireisl, Elsevier: Amsterdam, The Netherlands, 2005. 

[24] G.-Q. Chen, On nonlinear partial differential equations of mixed type, In: Proceedings of the 
4th International Congress of Chinese Mathematicians, Eds. L. Jin, K.-F. Liu, L. Yang, and 
S.-T. Yau, ISBN: 7-04-023267-7, Higher Education Press: Beijing, Jan. 2008. 

[25] G.-Q. Chen, J. Chen, and K. Song, Transonic nozzle flows and free boundary problems for 
the full Euler equations, J. Diff. Eqs. 229 (2006), 92-120. 

[26] G.-Q. Chen, J. Chen, and M. Feldman, Transonic shocks and free boundary problems for the 
full Euler equations in infinite nozzles, J. Math. Pures Appl. 88 (2007), 191-218. 

[27] G.-Q. Chen, C. M. Dafermos, M. Slemrod, and D. Wang, On two-dimensional sonic-subsonic 
flow, Commun. Math. Phys. 271 (2007), 635-647. 

[28] G.-Q. Chen and B.-X. Fang, Stability of transonic shock-fronts in steady potential flow past 
a perturbed cone, Discrete Conti. Dynamical Systems, 23, 85—114 (2009). 

[29] G.-Q. Chen and M. Feldman, Multidimensional transonic shocks and free boundary problems 
for nonlinear equations of mixed type, J. Amer. Math. Soc. 16 (2003), 461-494. 

[30] G.-Q. Chen and M. Feldman, Steady transonic shocks and free boundary problems in infinite 
cylinders for the Euler equations, Comm. Pure Appl. Math. 57 (2004), 310-356; 

[31] G.-Q. Chen and M. Feldman, Free boundary problems and transonic shocks for the Euler 
equations in unbounded domains, Ann. Scuola Norm. Sup. Pisa CI. Sci. (5), 3 (2004), 827- 
869. 

[32] G.-Q. Chen and M. Feldman, Existence and stability of multidimensional transonic flows 
through an infinite nozzle of arbitrary cross-sections, Arch. Rational Mech. Anal. 184 (2007), 
185-242 

[33] G.-Q. Chen and M. Feldman, Potential theory for shock reflection by a large-angle wedge, 

Proc. National Acad. Sci. USA (PNAS), 102 (2005), 15368-15372. 
[34] G.-Q. Chen and M. Feldman, Global solutions to shock reflection by large-angle wedges, Ann. 

Math. 2008, 106 pages (to appear; accepted on October 3, 2006). 
[35] G.-Q. Chen and M. Feldman, Regular shock reflection- diffraction and von Neumann sonic 

conjecture, In preparation, 2008. 
[36] G.-Q. Chen, M. Slemrod, and D. Wang, Vanishing viscosity method for transonic flow, Arch. 

Rational Mech. Anal. 189 (2008), 159-188. 
[37] G.-Q. Chen and H. Yuan, Uniqueness of transonic shock solutions in a duct for steady 

potential flow, J. Diff. Eqs. 2008 (to appear). 
[38] G.-Q. Chen, Y. Zhang, and D. Zhu, Existence and stability of supersonic Euler flows past 

Lipschitz wedges, Arch. Rational Mech. Anal. 181 (2006), 261-310. 
[39] S.-X. Chen, Linear approximation of shock reflection at a wedge with large angle, Commun. 

Partial Diff. Eqs. 21 (1996), 1103-1118. 
[40] S.-X. Chen, Asymptotic behavior of supersonic flow past a convex combined wedge, Chinese 

Annals Math. 19B (1998), 255-264. 
[41] S.-X. Chen, Supersonic flow past a concave double wedge, Sci. China, 41A (1998), 39-47. 
[42] S.-X. Chen, Global existence of supersonic flow past a curved convex wedge, J. Partial Diff. 

Eqs. 11 (1998), 73-82. 

[43] S.-X. Chen, Stability of transonic shock fronts in two-dimensional Euler systems, Trans. 

Amer. Math. Soc, 357 (2005), 287-308. 
[44] S.-X. Chen, Existence of stationary supersonic flows past a point body, Arch. Rational Mech. 

Anal. 156 (2001), 141-181. 
[45] S.-X. Chen, Stability of a Mach configuration, Comm. Pure Appl. Math. 59 (2006), 1-35. 
[46] S.-X. Chen, Mach configuration in pseudo-stationary compressible flow, J. Amer. Math. Soc. 

21 (2008), 63-100. 

[47] S.-X. Chen and B. Fang, Stability of transonic shocks in supersonic flow past a wedge, J. 
Diff. Eqs. 233 (2007), 105-135. 



SHOCK REFLECTION-DIFFRACTION 



23 



[48] S.-X. Chen, Z. Xin, and H. Yin, Global shock waves for the supersonic flow past a perturbed 
cone, Commun. Math. Phys. 228 (2002), 47-84. 

[49] Chen, S. and Yuan, H., Transonic shocks in compressible flow passing a duct for three- 
dimensional Euler systems, Arch. Rational Mech. Anal. 187 (2008), 523-556. 

[50] A. J. Chorin, Vorticity and Turbulence, Springer- Verlag: New York, 1994. 

[51] J. D. Cole and L. P. Cook, Transonic Aerodynamics, North-Holland: Amsterdam, 1986. 

[52] P. Colella and H. M. Glaz, Efficient solution algorithms for the Riemann problem for real 
gases, J. Comp. Phys. 59 (1985), 264-289. 

[53] P. Colella and L. F. Henderson, The von Neumann paradox for the diffraction of weak shock 
waves, J. Fluid Mech. 213 (1990), 71-94. 

[54] R. Courant, and K. O. Friedrichs, Supersonic Flow and Shock Waves, Springer- Verlag: New 
York, 1948. 

[55] C. M. Dafcrmos, Hyperbolic Conservation Laws in Continuum Physics, 2nd Ed., Springcr- 
Verlag: Berlin. 

[56] P. Daskalopoulos and R. Hamilton, The free boundary in the Gauss curvature flow with flat 
sides, J. Reine Angew. Math. 510 (1999), 187-227. 

[57] A. Yu. Dem'yanov and A. V. Panasenko, Numerical solution to the problem of the diffraction 
of a plane shock wave by a convex corner, Fluid Dynamics, 16 (1981), 720-725 (Translated 
from the original Russian). 

[58] R L. Deschambault and I. I. Glass, An update on non- stationary oblique shcok-wave reflec- 
tions: actual isopyenics and numerical experiments, J. Fluid Mech. 131 (1983), 27—57. 

[59] G. Dong, Nonlinear Partial Differential Equations of Second Order, Transl. Math. Mono- 
graphs, 95, AMS: Providence, RI, 1991. 

[60] V. Elling, Counterexamples to the sonic criterion, Preprint 2008. 

[61] V. Elling and T.-P. Liu, Supersonic flow onto a solid wedge, Comm. Pure Appl. Math. 2008 
(to appear). 

[62] V. Elling and T.-P. Liu, The ellipticity principle for self-similar potential flows, J. Hyper. 

Diff. Eqs. 2 (2005), 909-917. 
[63] B.-X. Fang, Stability of transonic shocks for the full Euler system in supersonic flow past a 

wedge, Math. Methods Appl. Sci. 29 (2006), 1-26. 
[64] R. Finn and D. Gilbarg, Three-dimensional subsonic flows, and asymptotic estimates for 

elliptic partial differential equations, Acta Math. 98 (1957), 265—296. 
[65] J. Francheteau and G. Mctivier, Existence de chocs faibles pour des systemes quasi-lineaires 

hyperboliques multidimensionnels, Asterisquc, 268, SMF: Paris, 2000. 
[66] I. M. Gamba and C. S. Morawetz, A viscous approximation for a 2-D steady semiconductor 

or transonic gas dynamic flow: existence theorem for potential flow, Comm. Pure Appl. Math. 

49 (1996), 999-1049. 

[67] I. M. Gamba, R. Rosales, and E. Tabak, Constraints on possible singularities for the unsteady 

transonic small disturbance (UTSD) equations, Comm. Pure Appl. Math. 52 (1999), 763—779. 
[68] D. Gilbarg and N. Trudinger, Elliptic Partial Differential Equations of Second Order, 2nd 

Ed., Springer- Verlag: Berlin, 1983. 
[69] I. I. Glass, Some aspects of shock-wave research, AIAA J. 25 (1987), 214-229. 
[70] I. I. Glass, Shock Waves and Man, The University of Toronto Press: Toronto, 1974. 
[71] H. M. Glaz, P. Colella, I. I. Glass, and R. L. Deschambault, A numerical study of oblique 

shock-wave reflection with experimental comparisons, Proc. Roy. Soc. Lond. A398 (1985), 

117-140. 

[72] H. M. Glaz, P. Colella, I. I. Glass, and R. L. Deschambault, A detailed numerical, graphi- 
cal, and experimental study of oblique shcok wave reflection, Lawrence Berkeley Laboratory 
Report LBL-20033 (1985). 

[73] H. M. Glaz, P. A. Walter, I. I. Glass, and T. C. J. Hu, Oblique shock wave reflections in SFq: 
A comparison of calculation and experiment, AIAA J. Porg. Astr. Aero. 106 (1986), 359—387. 

[74] J. Glimm, Nonlinear waves: overview and problems, In: Multidimensional Hyperbolic Prob- 
lems and Computations (Minneapolis, MN, 1989), 89-106, IMA Vol. Math. Appl. 29, Springer: 
New York, 1991. 

[75] J. Glimm, X.-M. Ji, J.-Q. Li, X.-L. Li, P. Zhang, T. Zhang, and Y.-X. Zheng, Transonic 
shock formation in a rarefaction Riemann problem for the 2-D compressible Euler equations, 
SIAM J. Appl. Math. 2008 (in press). 



24 



GUI-QIANG CHEN AND MIKHAIL FELDMAN 



[76] J. Glimm, C. Klingenberg, O. McBryan, B. Plohr, D. Sharp, and S. Yaniv, Front tracking 

and two-dimensional Riemann problems, Adv. Appl. Math. 6 (1985), 259—290. 
[77] J. Glimm and A. Majda, Multidimensional Hyperbolic Problems and Computations, Springcr- 

Verlag: New York, 1991. 
[78] J. Glimm, H. C. Kranzer, D. Tan, and F. M. Tangerman, Wave fronts for Hamilton- J acobi 

equations: the general theory for Riemann solutions in R n , Commun. Math. Phys. 187 (1997), 

647-677. 

[79] C.-H. Gu, A method for solving the supersonic flow past a curved wedge (in Chinese), Fudan 
Univ. J. 7 (1962), 11-14. 

[80] K. G. Gudcrley, The Theory of Transonic Flow, Pergamon Press: Oxford-London-Paris- 
Frankfurt; Addison- Wesley Publishing Co. Inc.: Reading, Mass. 1962. 

[81] O. Gues, G. Metivier, M. Williams, and K. Zumbrun, Existence and stability of multidimen- 
sional shock fronts in the vanishing viscosity limit, Arch. Rational Mech. Anal. 175 (2005), 
151-244. 

[82] O. Gues, G. Metivier, M. Williams, and K. Zumbrun, Navier-Stokes regularization of multi- 
dimensional Euler shocks, Ann. Sci. Ecole Norm. Sup. (4) 39 (2006), 75-175. 

[83] Jacques Hadamard, Legons sur la Propagation des Ondes et les Equations de 
I'Hydrodynamique, Hermann, Paris 1903 (Reprinted by Chelsea 1949). 

[84] E. Harabctian, Diffraction of a weak shock by a wedge, Comm. Pure Appl. Math. 40 (1987), 
849-863. 

[85] R. G. Hindman, P. Kutler, and D. Anderson, A two-dimensional unsteady Euler- equation 
solver for flow regions with arbitrary boundaries, AIAA Report, 79-1465 (1979). 

[86] H. G. Hornung, H. Oertel, and R. J. Sandeman, Transition to Mach reflexion of shock waves 
in steady and pseudodteady flow with and withiut relaxtion, J. Fluid mech. 90 (1979), 541-560. 

[87] J. K. Hunter, Transverse diffraction of nonlinear waves and singular rays, SIAM J. Appl. 
Math. 48 (1988), 1-37. 

[88] J. Hunter, Hyperbolic waves and nonlinear geometrical acoustics, In: Transactions of the 6th 
Army Conference on Appl. Math, and Computing, Boulder CO (1989), 527-569. 

[89] J. Hunter, Nonlinear wave diffraction, In: Geometrical Optics and Related Topics (Cortona, 
1996), 221-243, Progr. Nonlinear Diff. Eqs. Appl. 32, Birkhauscr Boston: Boston, MA, 1997. 

[90] J. K. Hunter and J. B. Keller, Weak shock diffraction, Wave Motion, 6 (1984), 79-89. 

[91] J. K. Hunter and A. Tesdall, Self-similar solutions for weak shock reflection, SIAM J. Appl. 
Math. 63 (2002), 42-61. 

[92] M. S. Ivanov, D. Vandromme, V. M. Formin, A. N. Kudryavtsev, A. Hadjadj, and D. V. 
Khotyanovsky, Transition between regular and Mach reflection of shock waves: new numerical 
and experimental results, Shock Waves, 11 (2001), 199-207. 

[93] J. B. Keller and A. A. Blank, Diffraction and reflection of pulses by wedges and corners, 
Comm. Pure Appl. Math. 4 (1951), 75-94. 

[94] D. Kindcrlchrer and L. Nirenberg, Regularity in free boundary problems, Ann. Scuola Norm. 
Sup. Pisa CI. Sci. (4), 4 (1977), 373-391. 

[95] V. P. Korobcinikov, Ed. Nonstationary Interactions of Shock and Detonation Wave in Gases, 
Nauka, Moscow, USSR, 1986 (in Russian). 

[96] V. P. Krehl and M. van der Geest, The discovery of the Mach reflection effect and its demon- 
stration in an auditorium, Shock Wave, 1 (1991), 3-15. 

[97] A. Kurganov and E. Tadmor, Solution of two-dimensional Riemann problems for gas dynam- 
ics without Riemann problem solvers, Numer. Methods Partial Diff. Eqs. 18 (2002), 584-608. 

[98] P. Kutler and V. Shankar, Diffraction of a shock wave by a compression corner, Part I: 
regular reflection, AIAA J. 15 (1977), 197-203. 

[99] P. D. Lax and X.-D. Liu, Solution of two-dimensional Riemann problems of gas dynamics 
by positive schemes, SIAM J. Sci. Comput. 19 (1998), 319-340. 

[100] R. J. LcVequc, Numerical Methods for Conservation Laws, Second Ed., Birkhauser-Verlag: 
Basel, 1992. 

[101] J. Li, T. Zhang, and S. Yang, The Two- Dimensional Riemann Problem in Gas Dynamics, 

Longman (Pitman Monographs 98): Essex, 1998. 
[102] T.-T. Li, On a free boundary problem, Chinese Ann. Math. 1 (1980), 351-358. 
[103] W.-C. Lien and T.-P. Liu, Nonlinear stability of a self-similar 3-dimensional gas flow, 

Commun. Math. Phys. 204 (1999), 525-549. 



SHOCK REFLECTION-DIFFRACTION 



25 



104] M J. Lighthill, The diffraction of a blast, I: Proc. Royal Soc. London, 198A (1949), 454-470; 

II: Proc. Royal Soc. London, 200A (1950), 554-565. 
105] F. H. Lin and L.-H. Wang, A class of fully nonlinear elliptic equations with singularity at 

the boundary, J. Gcom. Anal. 8 (1998), 583-598. 
106] T.-P. Liu, Hyperbolic and Viscous Conservation Laws, CBMS-NSF Regional Conf. Series in 

Appl. Math. 72, SIAM: Philadelphia, 2000. 
107] G. D. Lock and J. M. Dewey, An experimental investigation of the sonic criterion for 

transition from regular to Mach reflection of weak shock waves, Experiments in Fluids, 7 

(1989), 289-292. 

108] V. N. Lyakhov, V. V. Podlubny, and V. V. Titarcnko, Influence of Shock Waves and Jets 

on Elements of Structures, Mashinostroenie: Moscow, 1989 (in Russian). 
109] E. Mach, Uber den verlauf von Junkenwellenin der ebene und im raume, Sitzungsber. Akad. 

Wiss. Wicn, 78 (1878), 819-838. 
110] A. Majda and E. Thomann, Multidimensional shock fronts for second order wave equations, 

Comm. Partial Diflf. Eqs. 12 (1987), 777-828. 
Ill] A. Majda, Compressible Fluid Flow and Systems of Conservation Laws in Several Space 

Variables, Springer- Verlag: New York, 1984. 
112] A. Majda, One perspective on open problems in multi- dimensional conservation laws, In: 

Multidimensional Hyperbolic Problems and Computations (Minneapolis, MN, 1989), 217—238, 

IMA Vol. Math. Appl. 29, Springer: New York, 1991. 
113] A. Majda and A. Bertozzi, Vorticity and Incompressible Flow, Cambridge University Press: 

Cambridge, 2002. 

114] G. Metivier, Stability of multi- dimensional weak shocks, Comm. Partial Diff. Eqs. 15 (1990), 
983-1028. 

115] C. S. Morawetz, On a weak solution for a transonic flow problem, Comm. Pure Appl. Math. 
38 (1985), 797-818. 

116] C. S. Morawetz, Potential theory for regular and Mach reflection of a shock at a wedge, 

Comm. Pure Appl. Math. 47 (1994), 593-624. 
117] C. S. Morawetz, On steady transonic flow by compensated compactness, Meth. Appl. Anal. 

2 (1995), 257-268. 

118] H. Rcichcnbach, Contribtuon of Ernst Mach to fluid mechanics, Ann. Rev. Fluid Mech. 15 
(1983), 1-28. 

119] D. G. Schacflcr, Supersonic flow past a nearly straight wedge, Duke Math. J. 43 (1976), 
637-670. 

120] G. P. Schncycr, Numerical simulation of regular and Mach reflection, Phy. Fluids, 18 (1975), 
1119-1124. 

121] C. W. Schulz-Rinne, J. P. Collins, and H. M. Glaz, Numerical solution of the Riemann 

problem for two-dimensional gas dynamics, SIAM J. Sci. Comput. 14 (1993), 1394—1414. 
122] D. Serre, Perfect fluid flow in two independent space variables: Reflection of a planar shock 

by a compressive wedge, Arch. Rational Mech. Anal. 132 (1995), 15-36. 
123] D. Serre, Shock reflection in gas dynamics, In: Handbook of Mathematical Fluid Dynamics, 

Vol. 4, pp. 39-122, Eds: S. Friedlander and D. Serre, Elsevier: North-Holland, 2007. 
124] D. Serre, Spiral waves for the two-dimensional Riemann problem for a compressible fluid, 

Ann. Fac. Sci. Toulouse Math. (6), 5 (1996), 125-135. 
125] D. Serre, Multidimensional shock interaction for a Chaplygin gas, Arch. Rational Mech. 

Anal. 2008 (to appear), 10.1007/s00205-008-0110-z (Online). 
126] V. Shankar, P. Kutler, and D. Anderson, Diffraction of a shock wave by a compression 

corner, Part II: single Mach reflection, AIAA J. 16 (1978), 4—5. 
127] W.-C. Shcng and G. Yin, Transonic shock and supersonic shock in the regular reflection of 

a planar shock, Z. Angcw. Math. Phys. (ZAMP), DOI 10.1007/s00033-008-8003-4 (Online). 
128] M. Shiftman, On the existence of subsonic flows of a compressible fluid, J. Rational Mech. 

Anal. 1 (1952), 605-652. 
129] B. Skews and J. Ashworth, The physical nature of weak shock wave reflection, J. Fluid 

Mech. 542 (2005), 105-114. 
130] K. Song, The pressure- gradient system on non-smooth domains, Commun. Partial. Diff. 

Eqs. 28 (2003), 199-221. 
131] R. S. Srivastava, Starting point of curvature for reflected diffracted shock wave, AIAA J. 33 

(1995), 2230-2231. 



26 



GUI-QIANG CHEN AND MIKHAIL FELDMAN 



[132] R. S. Srivastava, Effect of yaw on the starting point of curvature for reflected diffracted 
shock wave (subsonic and sonic cases), Shock Wave, 17 (2007), 209-212. 

[133] B. Sturtevant and V. A. Kulkarny, The focusing of weak shock waves, J. Fluid. Mech. 73 
(1976), 1086-1118. 

[134] E. Tabak and R. Rosales, Focusing of weak shock waves and the von Neumann paradox of 

oblique shock reflection, Phys. Fluids, 6 (1994), 1874-1892. 
[135] A. M. Tesdall and J. K. Hunter, Self-similar solutions for weak shock reflection, SIAM J. 

Appl. Math. 63 (2002), 42-61. 
[136] A. M. Tesdall, R. Sanders, and B. Keyfitz, The triple point paradox for the nonlinear wave 

system, SIAM J. Appl. Math. 67 (2006), 321-336. 
[137] A. M. Tesdall, R. Sanders, and B. Keyfitz, Self-similar solutions for the triple point paradox 

in gasdynamics, SIAM J. Appl. Math. 2008 (to appear). 
[138] R. Timman, Unsteady motion in transonic flow, In: Symposium Transsonicum (IUTAM, 

Aachen, Sept. 1962), pp. 394-401, cd. K. Oswatitsch, Springer- Verlag: Berlin, 1964. 
[139] M. Van Dyke, An Album of Fluid Motion, The Parabolic Press: Stanford, 1982. 
[140] E. Vasil'cv and A. Kraiko, Numerical simulation of weah shock diffraction over a wedge 

under the von Neumann paradox conditions, Comput. Math. Math. Phys. 39 (1999), 1335- 

1345. 

[141] T. von Karman, The similarity law of transonic flow, J. Math. Phys. 26 (1947), 182-190. 
[142] J. von Neumann, Oblique reflection of shcoks, Navy Department, Bureau of Ordance, Ex- 
plosive Research Report, No. 12, 1943. 
[143] J. von Neumann, Collect Works, Vol. 5, Pergamon: New York, 1963. 

[144] P. Woodward and P. Colella, The numerical simulation of two-dimensional fluid flow with 

strong shocks, J. Comp. Phys. 54 (1984), 115-173. 
[145] P. Woodward, Simulation of the Kelvin- Helmholtz instability of a supersonic slip surface 

with the piecewise-parabolic method (PPM), In: Numerical Methods for the Euler Equations 

of Fluid Dynamics (Rocqucncourt, 1983), 493-508, SIAM: Philadelphia, PA, 1985. 
[146] X.-M. Wu, Equations of Mathematical Physics, Higher Education Press: Beijing, 1956 (in 

Chinese). 

[147] C. Xie and Z. Xin, Global subsonic and subsonic-sonic flows through infinitely long nozzles, 

Indiana Univ. Math. J. 56 (2007), 2991-3023. 
[148] Z. Xin and H. Yin, Transonic shock in a nozzle: two-dimensional case, Comm. Pure Appl. 

Math. 58 (2005), 999-1050. 
[149] G.-J. Yang, The Euler- Poisson-Darboux Equations, Yunnan University Press: Yunnan, 1989 

(in Chinese). 

[150] H. Yuan, On transonic shocks in two-dimensional variable-area ducts for steady Euler sys- 
tem, SIAM J. Math. Anal. 38 (2006), 1343-1370. 
[151] H. Yuan, Transonic shocks for steady Euler flows with cylindrical symmetry, Nonlinear 

Analysis, 66 (2007), 1853-1878. 
[152] G. I. Zabolotskaya and R. V. Khokhlov, Quasi-plane waves in the nonlinear acoustics of 

confined beams, Sov. Phys. -Acoustics, 15 (1969), 35—40. 
[153] A. R. Zakharian, M. Brio, J. K. Hunter, and G. M. Webb, The von Neumann paradox in 

weak shock reflection, J. Fluid Mech. 422 (2000), 193-205. 
[154] T. Zhang and Y. Zheng, Conjecture on the structure of solutions of the Riemann problem 

for two-dimensional gas dynamics, SIAM J. Math. Anal. 21 (1990), 593-630. 
[155] Y.-Q. Zhang, Steady supersonic flow past an almost straight wedge with large vertex angle, 

J. Diff. Eqs. 192 (2003), 1-46. 
[156] Y. Zheng, Systems of Conservation Laws: Two- Dimensional Riemann Problems, 

Birkhauscr: Boston, 2001. 
[157] Y. Zheng, A global solution to a two-dimensional Riemann problem involving shocks as free 

boundaries, Acta Math. Appl. Sin. Engl. Ser. 19 (2003), 559-572. 
[158] Y. Zheng, Two-dimensional regular shock reflection for the pressure gradient system of 

conservation laws, Acta Math. Appl. Sin. Engl. Ser. 22 (2006), 177-210. 



SHOCK REFLECTION-DIFFRACTION 27 
Department of Mathematics, Northwestern University, Evanston, IL 60208-2730, 

USA 

E-mail address: gqchenSmath.northwestern.edu 

Department of Mathematics, University of Wisconsin, Madison, WI 53706-1388, USA 
E-mail address: feldmanamath.wisc.edu 



